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Abstract 

Surface effects are generally prevailing in confined colloidal systems. Here we report on dispersed 
nanoparticles close to a fluid membrane. Exact results regarding the static organization are derived 
for a dilute solution of non-adhesive colloids. It is shown that thermal fluctuations of the membrane 
broaden the density profile, but on average colloids are neither accumulated nor depleted near the 
surface. The radial correlation function is also evaluated, from which we obtain the effective pair- 
potential between colloids. This entropically-driven interaction shares many similarities with the 
familiar depletion interaction. It is shown to be always attractive with range controlled by the 
membrane correlation length. The depth of the potential well is comparable to the thermal energy, 
but depends only indirectly upon membrane rigidity. Consequenses for stability of the suspension 
are also discussed. 

PACS numbers: 05.20.-y, 82.70.Dd, 87.16.Dg 
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I. INTRODUCTION 

Fluid membranes are soft surfaces, self-assembled from surfactant solutions They 
assume a large variety of shapes and topologies, which have been accurately explained in 
terms of bending energy 0, |J. In most practical realizations, however, membrane suspen- 
sions are not pure but incorporate colloidal entities as well. In living systems for instance, 
lipid bilayers organize the cell into compartments that keep apart different chemical envi- 
ronments [6]. Biological membranes are therefore in contact with various kinds of proteins, 
macro-ions or more complex structures ^. Membrane phases used in detergency or cos- 
metics also include numerous additives, like macromolecules or colloids, in order to improve 
efficiency and to control viscoelastic properties j^. One is then naturally inclined to inves- 
tigate transformations that may occur upon addition of colloids, and several studies have 
recently been devoted to these complex systems. Essentially, these have focused on the 

n 

softening of membrane resulting from the depletion of spherical and rod-like colloids 7], or 



the depletion 



m 



and adsorption 
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, |l^ of flexible polymers, but the predicted effects 
are usually small compared to the bare rigidity of the bilayer. Bending of a membrane upon 
colloid adsorption has also drawn growing interest because of potential applications for drug 



encapsulation and gene delivery 



L growing interest b( 



Nevertheless, a common feature shared 
by most theoretical studies is that membrane fluctuations are systematically disregarded. 
This point clearly illustrates the technical difficulty to couple bulk and surface degrees of 
freedom. 

Generally, the mutual influence of bulk and surface properties on each other is a chal- 

The situation has now been clarifled for bidispersed hard-sphere 
|21|, |2^: when in contact with a flat substrate, excluded- volume effects are 
known to push the larger beads toward the wall of the sample Q]. Recent experiments 
done with curved or corrugated surfaces have shown that geometric features of the surface 
ca. a.o ceate and .nodulate ...epic fo.ce «e.d. Q Q. These depletion feces can be 
used to grow oriented colloidal crystal, with numerous potential applications such as the 
fabrication of photonic bandgap crystals [2^. The theoretical description of those systems 
usually requires advanced density functional techniques |2l|, that have been adapted to 



lenging problem 118|, Ui 
2],! 



suspensions 



study depletion potentials close to arbitrarily shaped substrates 



27| . In this work, we follow 



a different line and present some new flndings regarding the static organization of nanopar- 
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tides near a fluctuating surface. Given the increasing complexity of the problem, we focus 
on the simplest system consisting of a monodispersed, dilute solution of non-adhesive col- 
loids. This allows us to derive exact results and to highlight non-trivial phenomena, such as 
membrane-induced interaction between the colloids. 

This paper is organized as follows. In Sec. IHl we present an exact computation of the 
partition function of the global system. The determination of the variation of the particle 
density upon distance is the aim of Sec. IIIII The computations of the radial distribution 
function and the resulting pair-potential and related discussion are presented in Sec. IIVI 
We draw some concluding remarks in the last section. Finally, some technical details are 
relegated to the appendices. 

II. THE PARTITION FUNCTION 

Consider a fluctuating membrane in contact with a colloidal suspension, consisting of 
nanoparticles immersed in some solvent. We assume that the colloids cannot permeate 
through the bilayer on experimental time-scales, in such way that the membrane act as 
a flexible but impenetrable wall for the particles. In this paper, we shall use the notation 
r = (p, z), where p = (x, y) is the transverse vector and z the perpendicular component. The 
position of the (almost flat) membrane is specified through the displacement field h{x,y). 
The surface fluctuates around the horizontal plane z = 0, so that the height h [x, y) may 
take either positive or negative values. The equilibrium statistical mechanics of membranes 
is based on the Helfrich Hamiltonian l2l 



where {K,a,fi) are the elastic constants of the membrane. In what follows, we shall rather 
use the rescahng parameters: k, = (3k, a = [3a, and /i = where we define as usually 
(3 = l/ZcfiT, with T the absolute temperature and ks the Boltzmann constant. Let us also 
introduce the height-height correlation function 



from which one obtains the mean-squared fluctuations = G{0). Some properties of G 
are recalled in Appendix 1X1 




(1) 



G (p - pO = (Mp) ))o - (Mp))o (MpO) 



' 



(2) 
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We intend to compute some statistical properties of the system under investigation, 

1. For a dilute so- 



namely the particle density profile and radial distribution function 
lution, one may ignore the interactions between particles and treat them as an ideal gas 
in confined geometry. Since we are interested in those systems containing a considerable 
number of particles (thermodynamic limit), the physical quantities are independent of the 
particular choice of statistical ensemble. For convenience, we shall consider both the canon- 
ical and the grand canonical ensembles. The grand canonical partition function, denoted 

Zqi is the Laplace transform of the canonical partition function, that is 

+00 

Za = Y. {N) , (3) 

N=0 

with / the fugacity. In the above equality, the canonical partition function Zc (N) is 

Z. (N) = I dh,...dhr. I ^(-) . (4) 

The functional integral extends over all configurations of the field h, weighted with the 
Helfrich energy (P). The positions Rat = (ri, ...,vn) of the particles are restricted to the 
upper side of the three-dimensional space limited by the membrane. In the above equation, 
the thermal wavelength A results from integration over the particles momenta, and colloid- 
membrane interactions are accounted for through the (contact) hard core potential V (r) 

, , ii Zi> h (xi, Vi) , 
-00 , otherwise . 

With these considerations, we first examine the canonical partition function. It can be 
written 

- f f[ dh.ipN (Rn) , (6) 
■ ^ i=i 



^ . N 

z = 

where it is convenient to define the function (p^ 



r ^ 

ifr, (Riv) = Zo' X / Vhe-^'^^''^ l[e [z, - h (x„ y,)] . (7) 

i=i 

Here, Zq = J Vhexp {—/STi. [h]} is the partition function of a membrane in the absence of 
particles (A^ = 0), and 9 (x) is the step function. It is easy to see, from its definition 
that ifN satisfies the following boundary conditions at infinity 

y^^N {pi, ■■; Pn; zi, zn)\,,=,„=,^=^oo = > (8) 
ifN {pi, Pn; zi, zn)\,^^,,,^,^^+^ = 1 , (9) 



for fixed values of the transverse vectors (pi, p^v). The function at (Rat) has to be 
understood as the partition function of a membrane whose configurations are subjected to 
N restrictions h{xi,yi) < Zi, with 1 < i < A^. It is generally a complicated function of 
the relative transverse distances \pi — pj\ and perpendicular components Zi. Indeed, the 
translation symmetry is preserved in the parallel directions, but due to the presence of the 
membrane, this symmetry is broken in the perpendicular one. After some algebra detailed 
in Appendix IHl we find that v^at is given by 



dz[... / dz'^exp -- J2 <■ [Gn% -4 ' (10) 

■oo J-oo [ ^ ^ -^^ J 



where ^at is the squared matrix of order A^ whose elements are the propagators [QN]ij = 
G {pi — pj) defined in Eq. 0. This expression is compatible with the boundary conditions 
(jHl) and Q, and more generally we can deduce from Eq. (jl(J|) the fundamental property 
according to which 

< V^Af (Pl, PN] Zl, Zat) < 1 , (11) 

for all values of the position vectors (ri, fat). The right-hand side inequality relies on the 
fact that y^Af is an incomplete Gaussian multiple-integral. 

Now, we direct our attention to the grand canonical function. Formally, it writes as 

+00 ^ r r ^ 

= E/'^A^ / ^f'^''"''"' / lld'-^0[z,^-hip,)] 

N=Q ' -^1=1 

= Jvhe-'^'^^''^exp!^-^jdhe[z-h{p)]^. (12) 
Integrating over the ^— variable then yields 

Zg = Zoe^^/2A ^ J _ (^3) 

^0 

Here, Q = S xL is the total volume occupied by the system under investigation, S is the area 
of the horizontal plane and L/2 represents the upper bound of the perpendicular coordinate 
z. The above functional integral can be easily calculated, and we find 

= exp {fSG (0) /2A^} , (14) 
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with G (0) = l/yU the Fourier transform of the propagator G at g = 0. To obtain this 
result, we have apphed the general formula ()B5|) of Appendix B to the particular source 
J (p) = —fX^^. Finally, we have the simple, exact expression for the partition function Zq 

2c = 2.xexp{g}xexp{0}. (15) 

Noticeably, splits into three parts. The first two factors are the standard contributions: 
Zo is the partition function of the membrane in a particle-free solvent, and exp {fQ/2X^} is 
the partition function of an ideal gas of colloids (remember that the particles are restricted 
to the upper half of the space). All the information concerning the interplay between bulk 
and surface contributions finally factorizes in the last term of Eq. ()15|) . 

With this expression of Zq, we are now able to evaluate the average number of colloids 

In addition to the usual term fQ/2X^, we find a second contribution that happens to 
be negligible at low concentration (/<<!) or for very large perpendicular extension 
[L >> /i~^A~'^). For the sake of simplicity, we will assume thereafter that these require- 
ments are fulfilled so that the bulk concentration is 

P»^§-/A-. (17) 

Note, however, that the second contribution in Eq. (fTHIl could not be disregarded in a strongly 
confining system. This point might be relevant for experimental realizations involving col- 
loids in a lamellar phase or in a sponge phase 31 1 of membranes. 



III. PARTICLE DENSITY PROFILE 

The concentration profile of an ideal gas of colloids in contact with a rigid wall located at 
2; = is simply pnw = Poo x ^(-2)- For a fiexible interface, the situation is quite different as 
thermal undulations are expected to broaden the distribution. In this section, we examine 
the mean-value of the particle density at point r defined by 

p(r) = /^53(r-r,)\=Ar(53(r-rO) , (18) 
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with ^3 the three-dimensional Dirac distribution. The averaging procedure imphes integra- 
tion over the colloid configurations as well as internal degrees of freedom of the membrane. 
Explicitly, we have 

+00 fNf c 



2, 



G 



The last functional integral can be computed making a simple translation of the field h — > 
h + //A^/i. We finally get 

p{r) = fX-'x^,{p;z + f/X'jl) , (20) 
where ipi is a particular function of type (jTU)), that is 

ip,{p-z + fl\^p) = [2TiGm'"'j dz'e^^l--G-\Q)z^\ . (21) 

For = 1, the squared matrix Qi in relation (jiup reduces to G (0) = Recalling that 
/A^^ equals the bulk density poo, we find the final expression for the density profile 

1 



P {z) = Poo X 
with the characteristic length 



(22) 



zq = PooP ^ • (23) 

The concentration profile is shown in fig. (0). As expected, it only depends on the perpendic- 
ular distance z (homogeneity property in the parallel directions). For fixed parameter p, the 
scale Zq becomes smaller as the particle density is decreased. The physical meaning of this 
length can be understood as follows. When in contact with the colloidal solution, the mem- 
brane experiences the osmotic pressure of the particles. At low concentration, this pressure 
is proportional to the concentration of particles in contact with a fiat surface Posm = ksTpoo. 
Indeed, we see from Eq. (fT!?|) that integration over the colloid positions leads to an effective 
Hamiltonian for the membrane Tie// = 'Ho + Posm J (P'phi^p). The average position of the 
membrane is then shifted to its new value {h) = —zq, so that the concentration profile (or 
the dividing surface) is translated from the same distance. For a symmetric system (particles 



on both sides with the same chemical potential), this length would just vanish. Note that 
when z = —zq, the local particle density is reduced to half of the bulk value, that is Poo/2. 

Now, let us compute the excess particle density, defined as the first moment of the density 
profile 3^ 

/ + CO 
dz [p {z) -pooX 0{z + z^)] . (24) 
-oo 

r is called the adsorption of the species, since a large positive value of F is an evidence 
for particle accumulation at the surface. Conversely, a negative value of F indicates that 
concentration in the surface vicinity is lower than concentration in the bulk phase. Using 
relation (j221), we find after a simple integration 

F = , (25) 

meaning that there is, on average, neither accumulation nor depletion of particles near 
the membrane. Actually, this result comes from the cancellation of two effects: insertion 
of particles into the holes and valleys of the rough surface exactly compensates for the 
depletion from the convex regions, as can be seen in fig. (^. Eq. ()25|) also implies that there 
is no additional contribution to the interfacial tension, and consequently no spontaneous 
curvature of the membrane induced by the colloids. However, this result is only valid for 
vanishing particle radius a. Although finite size effects are not easil y in cluded in the theory, 
one does not expect these results to hold for a ~ any longer [ifil Is^^. 

IV. RADIAL DISTRIBUTION FUNCTION AND EFFECTIVE POTENTIAL 

To better characterize the statistical properties of the system, we now focus on the pair 
distribution function defined as 

E,^,- (^3 (r - r*) ^3 (r'-r,)) 

g (r, rO = ^^^^ ^ ' ^ , ^ " (^6) 

P (r) P (r') 

Using the same techniques as before, we find that g (r, r') can be expressed in terms of the 
V?Ar functions defined in (ITUj) . Without further detail, we obtain the formal expression 

g (r r') = (p, p'; z + ^0, z' + 2:0) 

i^l{p]Z + Zo)ipi{p']Z' + Zq) ' 
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As one could expect, only yjj's with = 1 and 2 come out of the calculations. For the sake 
of completeness, we give explicitly the radial distribution function 



g (r, r 



(2vr)-Mdet ^,]-^/^ rf.i rf., exp {-1 (.1, ^,-1 (-) } 



l + erff^ 



l + erf(ilg) 



(28) 



where Q2 is the 2x2 correlation matrix 



G(0) G{\p-p'\) 
G{\p-p'\) G(0) 



(29) 



and det ^2 = G (0)' -G{\p- p'\y. Once again, the result depends on the relative transverse 
distance |p — whereas it varies with z and z' separately. We emphasize that this pair 
correlation function would be identically unity for an ideal gas confined by a rigid wall. 
Here, surface fluctuations give rise to hulk correlations between colloids over distances that 
depends both on their separation and on their upright distance from the membrane. For fixed 
z and z\ g (r, r') is maximum when p = p' and decreases to 1 as the transverse separation 
increases. The equality g (r, r') = 1 (no correlations) is only achieved either when z or z' 
goes to +00 at fixed parallel distance \p — or when \p — p'\ —>■ +00 at fixed z and z'. 
Indeed, the former requirement expresses that the colloids do not feel the surface anymore 
at elevations higher than whereas the latter asserts that correlations vanish at parallel 
separation much larger than the membrane correlation length (see Appendix A). 

With the help of the computed one- and two-points distribution functions, we can extract 
the membrane- induced interactions between particles. According to Eq. (jH)), an effective A^- 
body potential Uiji, .., r^r) may be defined through 



(PN {pl,..,pN]Zi, 



Zn) 



-/3W(ri,.,riv) 



(30) 



The analyze developped in this report indicates that the many-body interaction decomposes 

as 



W(ri, ...,Yn) = ^Wi(ri) + ^W2(ri,rj) + . . .+l(j^{ri, ..,rAr) . 



(31) 



Here, Ui is the effective external potential resulting from thermal undulations of the mem- 
brane. We easily find 



(r) = -kBTln 



Poo 



(32) 
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As the reduced density p{z) j p^o ranges from to 1, this potential is always repulsive 
and tends to move the particles away from the surface. Note that the corresponding force 
F = —dlAi/dz exerted on a particle by surface undulations remains finite at "contact": 
F{—zq) = ksT/^^ (recall that —zq is the average position of the membrane under the 
osmotic pressure of the colloids). Surprisingly, this force increases as the roughness of the 
surface decreases, but it has to be this way as F eventually diverges in the hard wall 
limit. 

Regarding the potential of mean force U2 (r, r'), we find 

U2{v,v') = -kBT\ng{\p- p'\,z,z') . (33) 

The normalization of g (r, r'), Eq. (jTTj) . ensures that only two-body terms are accounted for. 
The pair potential U2 is shown in fig. Q for a membrane with no surface tension (cr = 0) and 
for fixed z = z'\ it is negative at short parallel distances and vanishes at large separations. 
Accordingly, colloids that are close to the membrane tend to aggregate even if there are 
only hard-core repulsions in our description. We can evaluate the depth of this potential: 
diagonalizing the quadratic form in leads to 



U2 (r = r') = -ksT In 2 + keT In 



(34) 



Of course, the interaction still depends on the ^-position of the pair of colloids. The depth 
of the potential increases with decreasing altitude, and is of order ksT for z = —zq- At 
larger separation, U2 displays a tiny repulsive barrier (O.OlfcsT at most). In the limit 
d = \p — p'\ ^ , we show in Appendix C that the potential of mean force vanishes 
exponentially 

U2 (d) ~ kBT^e-'"^\\ , (35) 

for z = z' = Zq. 



V. DISCUSSION 



In this paper, we have adapted the usual many-body statistical mechanics in order to 
include an additional degree of freedom, namely the thermal fluctuations of the membrane. 
It has been shown that surface undulations broaden the density profile and generates corre- 
lations among an essentially ideal gas system. As can be seen in Eq. (fT^. partial integration 
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over the positions induces a linear coupling with the membrane height in the case of point- 
like particles. For finite-size objects, one would certainly expect further couplings with 
membrane curvature as well as with higher-order terms, but this point is far beyond the 
scope of this paper. 

As a consequence of surface fiuctuations, the colloids attract each other through the po- 
tential of mean force IA2- Interstingly, the only remaining signature of the elastic parameters 
of the membrane are the length scales ^y, and zq. For two particles at a given position, the 
configurational entropy of the membrane increases as the colloids come closer. This effective 
pair interaction is of order of the thermal energy for particles near the average position of 
the membrane, and is in many respects similar to the familiar depletion interaction. In 
particular, it would simply sum up with a direct colloid-colloid potential in a more realistic 
system. Note that our approach is to some extent peculiar, in the sense that we trace out 
the degrees of freedom of the "slow" variable, ending up with an effective Hamiltonian for 
the small particles. This procedure usually leads to a very poor description of the system, 
because one has generally to resort to uncontrolled approximations. Here, the situation is 
more refined as we managed to perform exact calculations. The potential of mean force 
is therefore expected to be very accurate for colloids much smaller than the correlation 
length . 

Finally, we would like to emphasize the fact that lAi and IA2 are the dominant interactions 
at low concentrations. The framework developed in this paper would in principle allow 
us to evaluate the relative weight of many-body contributions, but these would only be 
relevant at the onset of an hypothetical aggregation. Here however, the drift force from the 
membrane always move the particle away from the surface so that two-body attraction is 
not prevailing. One could still imagine to enforce particle accumulation near the membrane 
through a small attractive colloid-membrane interaction. Whether membrane fluctuations 
could then induce surface cristallization is an interesting point, work on this question is 
currently under progress. 
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APPENDIX A 

In this appendix, we recall a couple of results on fluid membranes. The height-height 
correlation function is defined as 

G{p) = (/i(p)MO))o-(Mp))o(MO))o 

'(^W+W^ ' ^^^^ 

with the notation q = |q|. Here, the thermal averages (. . .)q are performed with the Helfrich 
Hamiltonian (0) in the absence of particle. For a bilayer without surface tension (cr = 0), 
the integral over the Fourier modes leads to 

G (p) = --e^kei ( V2-^] , (A2) 



^ " V eii. 

with kei(a;) = Im [Kq (xe*'^'''^)] is a Kelvin function j^^. We then define the mean roughness 
of the membrane = (7(0)^^^ = 2~^/^(K./i)~^/^, and the in-plane correlation length = 
2^/^ (k//!)^'^^ characterizing the exponential decay of G (p) at large distances 

G (p) ~ e-'^/^ii , p»ei|. (A3) 

APPENDIX B 

The aim of this appendix is the proof of formula that defines the function (pN. To 
this end, we first compute its multiple derivative 



dzi...dzN 



Write the integral form of the Dirac distribution 



AT 



5 [z, - h (x„ y,)] = I _L^^A^.-^^^.v.)\ , (B2) 



oo 



allows us to rewrite Eq. (jBljl as 

oo j=l 



dzi...dzN 



27r 
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X ^ . (B3) 

^0 



Here, we have introduced the source 

N 

J[p) = -iY^kj52{p-p,) , (B4) 

with the two-dimensional vector pj = {xj,yj). The functional integration in relation ()B3jl is 
trivial, since the Hamiltonian 7i [h] is quadratic in the field h. We simply give the result 

I = ^^P[^J d'p J <fp'J (p) G{p- p') J ip') j , (B5) 

where G {p — p') is the membrane propagator, relation (j2I). Now, replace the source J (p) 
by its definition ()B4|) to find 

N 



\l d^pj d^p'J ip) G{p- p') J ip') = -i 5^ A;,G (p. - p,) k, 
J J j=i 



(B6) 



We therefore obtain 



dzi...dz 



/• + 00 ^ f 1 ^ 1 

= {27.)-''/'[detg^]-'/'xexpi^~^-f^z[.[g^\.z'^ . (B7) 

The N'^ coefficients [GN]ij = G {pi — Pj), I < i,j < N, define a squared matrix Qjy. Remark 
that the above multiple integral is Gaussian. A straightforward integration yields the explicit 
expression of the function 99 ^r, formula (fTU|). 



APPENDIX C 

At large transverse separation (i = |p — p'| ^ ^y, it is possible to evaluate completely the 
two-points function ip2- Indeed, the propagator satisfies in this limit G{d) <^ G{0) = so 
that we are naturally lead to define the small parameter a = G{d)/G{0). At first order, one 
has det^ = ^1(1 + 0{a^)), and 

(^1, ^2) h) =i + i- 2«7? + 0{a^) . (CI) 
Expanding the exponential in Eq. (fTUI) up to first order, we find 

/z+zo rz'+zo |' 2_|_ 2^/ \ 

dz, J dz^exp |-^^ j (1 + «^ + 0{a')j 

= v.i(z),.i(zO + ^exp{-iii^}exp{-i^^ (C2) 
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Taking the logarithm of this expression and using the definition ()A2|) for G{d) finally leads 
to Eq. dng). 
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Perpendicular distance 



FIG. 1: Reduced density p{z)/ poo as a function of the reduced perpendicular distance z/£^±. The 
density profile is symetric with respect of the dotted line z = —zq, so that the number of particles 
that enter the depressions of the undulating surface exactly compensate for the depleted ones. 
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Transverse distance 

FIG. 2: Equal-height, effective pair-potential pU2{p, p' , z = z') as a function of the reduced trans- 
verse distance \p — p'\ /.^|| . The curves corresponds to different values of the perpendicular distance. 
From bottom to top: {z + zo)/^± = 0, 0.5, and 1. The depth of the potential well is of order ksT 
for z + zq = 0, and vanishes exponentially at large separations. 
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